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ABSTRACT 

We address the dusty wind problem, from the point where dust formation has been 
completed and outward. Given grain properties, both radiative transfer and hydro- 
dynamics components of the problem are fully defined by four additional input pa- 
rameters. The wind radiative emission and the shape of its velocity profile are both 
independent of the actual magnitude of the velocity and are determined by just three 
dimensionless free parameters. Of the three, only one is always significant — for most 
of phase space the solution is described by a set of similarity functions of a single 
independent variable, which can be chosen as the overall optical depth at visual Ty. 
The self-similarity implies general scaling relations among mass loss rate (M), lumi- 
nosity (L) and terminal velocity (uoo)- Systems with different Af, L and v^a but the 
same combination M /L^^^ necessarily have also the same Mvoo/L. For optically thin 
winds we find the exact analytic solution, including the effects of radiation pressure, 
gravitation and (sub- and supersonic) dust drift. For optically thick winds we present 
numerical results that cover the entire relevant range of optical depths, and summa- 
rize all correlations among the three global parameters in terms of Ty. In all winds, 
M (xv^{l + tyY'^ with a proportionality constant that depends only on grain prop- 
erties. The optically thin end of this universal correlation, M oc u^, has been verified 
in observations; even though the wind is driven by radiation pressure, the luminosity 
does not enter because of the dominant role of dust drift in this regime. The M-L 
correlation is M oc (Lry)'^/^(1 -I- ry)°-^°^. At a fixed luminosity, M is not linearly 
proportional to Ty, again because of dust drift. The velocity- luminosity correlation 
is Vao oc (LTy)^/'*(l -f Ty)~°-*^^, explaining the narrow range of outflow velocities 

/ * 4/3 \ — 

displayed by dusty winds. Eliminating ry produces — AM {\ + B M / L\ , 

where A and B are coefficients that contain the only dependence of this universal 
correlation on chemical composition. At a given L, the maximal velocity of a dusty 
wind is Wmax oc L^/** attained at M oc L"^/^, with proportionality coefhcients derived 
from A and B. 

Key words: circumstellar matter — dust — infrared: stars — stars: AGB and post- 
AGB — stars: late-type — stars: winds, outflows 



1 INTRODUCTION 

Stars on the asymptotic giant branch (AGB) make a strong 
impact on the galactic environment. Stellar winds blown 
during this evolutionary phase are an important component 
of mass return into the interstellar medium and may account 
for a significant fraction of interstellar dust. These dusty 
winds reprocess the stellar radiation, shifting the spectral 
shape toward the infrared. There are good indications that 
they dominate the IR signature of normal elliptical galax- 
ies (Knapp, Gunn & Wynn- Williams, 1992). In addition to 



its obvious significance for the theory of stellar evolution, 
the study of AGB winds has important implications for the 
structure and evolution of galaxies. 

Because of the reddening, dusty winds are best studied 
in the infrared. Since the observed radiation has undergone 
significant processing in the surrounding dust shell, inter- 
pretation of the observations necessitates detailed radiative 
transfer calculations. The complexity of these calculations 
is compounded by the fact that the wind driving force is 
radiation pressure on the grains; complete calculations re- 
quire a solution of the coupled hydrodynamics and radiative 
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transfer problems. Traditionally these calculations involved 
a large number of input parameters that fall into two cat- 
egories. The first one involves the dust properties; widely 
employed quantities include the dust abundance, grains size 
distribution, solid density, condensation temperature and 
absorption and scattering efficiencies. The second category 
involves global properties, including the stellar temperature 
T, , luminosity L = 10''L4 L©, mass M = Mo Mq and the 
mass-loss rate M = 10~'^M-e Mq yr"\ 

The large number of input quantities complicates mod- 
eling efforts, making it unclear what are the truly indepen- 
dent parameters and which properties can actually be de- 
termined from a given set of data. In a previous study we 
noted that the dusty wind problem possesses general scal- 
ing properties such that, for a given type of grains, both 
the dynamics and radiative transfer depend chiefly on a sin- 
gle parameter — the overall optical depth (Ivezic & Elitzur 
1995, hereafter IE95). In a subsequent study we established 
in full rigor the scaling properties of the dust radiative trans- 
fer problem under the most general, arbitrary circumstances 
(Ivezic & Elitzur 1997, hereafter IE97). Here we extend rig- 
orous scaling analysis to the other aspect of the dusty wind 
problem, the dynamics. 

In scaling analysis the basic equations are reduced to 
the minimal number of free parameters that are truly inde- 
pendent of each other. By its nature, such analysis is driven 
by the underlying mathematics, and the proper free param- 
eters are not necessarily convenient for handling the data. 
Attempting to make our presentation tractable we have sep- 
arated it into two parts. The present paper discusses all 
the theoretical and mathematical aspects of the dusty wind 
problem. In a companion paper (Ivezic & Elitzur 2001, pa- 
per II hereafter) written as a stand-alone, the observational 
implications are discussed separately on their own. Read- 
ers mostly interested in practical applications may proceed 
directly to paper II. 



2 UNDERLYING THEORY 
2.1 Problem Overview 

The complete description of a dusty wind should start at its 
origin, the stellar atmosphere. Beginning with a full atmo- 
spheric model, it should incorporate the processes that initi- 
ate the outflow and set the value of M. These processes are 
yet to be identified with certainty, the most promising are 
stellar pulsation (e.g. Bowen 1989) and radiation pressure on 
the water molecules (e.g. Elitzur, Brown & Johnson 1989). 
Proper description of these processes should be followed by 
that for grain formation and growth, and subsequent wind 
dynamics. 

An ambitious program attempting to incorporate as 
many aspects of this formidable task as possible has been 
conducted over the past few years, yielding models in qual- 
itative agreement with observations (see Fleischer, Winters 
& Sedlmayr 1999 and Hofner 1999, and references therein). 
However, the complexity of this undertaking makes it diffi- 
cult to assess the meaning of its successes. In spite of con- 
tinuous progress, detailed understanding of atmospheric dy- 
namics and grain formation is still far from complete. When 
involved models succeed in spite of the many uncertain in- 



gredients they contain, it is not clear whether these ingredi- 
ents were properly accounted for or are simply irrelevant to 
the final outcome. 

Fortunately, the full problem splits naturally to two 
parts, as recognized long ago by Goldreich & Scoville (1976, 
GS hereafter). Once radiation pressure on the dust grains 
exceeds all other forces, the rapid acceleration to super- 
sonic velocities produces complete decoupling from the ear- 
lier phases that contain all the major uncertainties. Subse- 
quent stages of the outfiow are independent of the details 
of dust formation — they depend only on the final properties 
of the grains, not on how these grains were produced; the 
supersonic phase would be exactly the same in two differ- 
ent outflows if they have the same mass-loss rate and grain 
properties even if the grains were produced by entirely dif- 
ferent processes. Furthermore, these stages are controlled by 
processes that are much less dependent on detailed micro- 
physics, and are reasonably well understood. And since most 
observations probe only the supersonic phase, models de- 
voted exclusively to this stage should reproduce the observ- 
able results while avoiding the pitfalls and uncertainties of 
dust formation and the wind initiation. For these reasons, 
the GS approach with its focus on the supersonic phase has 
been widely used in studies of the dusty wind problem (in- 
cluding recent ones by Netzer & Elitzur 1993, NE hereafter, 
and Habing, Tignon & Tielens 1994, HTT hereafter). This 
is the problem we address here. 



2.1.1 Overall Plan 

We consider a spherical wind in steady state (the steady- 
state assumption is adequate as long as the wind structure 
is not resolved in too fine details; see IE95). Our starting 
point is the radius ri beyond which the properties of indi- 
vidual dust grains do not change and radiation pressure is 
the dominant force on the envelope. When positions in the 
shell are specified in terms of the scaled radius y = r/ri, 
the shell inner boundary is always at y = 1 and the actual 
magnitude of ri drops out of the problem. The equation of 
motion is pdv/dt — JT, where is the net outward radial 
force per unit volume, v is the gas velocity and p = nnmp 
is its density (nn is the number density of hydrogen nuclei 
and TUp is the proton mass). In steady-state dt = dr/v and 
the equation becomes 

^=2arr,, (1) 

where a^r — T I p is the acceleration associated with force T . 
Since arrx has dimensions of v , any force can always be 
characterized by the velocity scale it introduces. 

Winds of interest are highly supersonic, therefore the 
gas pressure gradient can be neglected (see also The 
expansion is driven by radiation pressure on the dust grains 
and is opposed by the gravitational pull of the star. The dust 
and gas particles are coupled by the internal drag force. For 
each of these three force components we first derive the char- 
acteristic velocity scale and the dimensionless profile associ- 
ated with its radial variation. With the resulting expressions 
we identify all the dimensionless free parameters and formu- 
late the problem in terms of independent dimensionless vari- 
ables, resulting in two coupled equations for the dynamics 
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Carbon 


Silicate 




0.1 


0.1 


Tc (K) 


800 


800 


Qv 


2.40 


1.15 


Q, 


.599 


.114 




5.97 


2.72 



Table 1. Standard parameters for dust grains used in all numer- 
ical calculations. The efficiency factors are from Hanner 1988 for 
amorphous carbon, and Ossenkopf, Henning & Mathis 1992 for 
(the "warm" version of) silicate grains. The grain size a and subli- 
mation temperature Tc are assumed. The lower part lists derived 
quantities; Qv is the efRciency factor for absorption at visual; Q, 
is the Planck average at the stellar temperature of the efficiency 
coefficient for radiation pressure (equation |4|); is defined in 
equation 



and radiative transfer. We proceed to solve the mathemat- 
ical problem, and afterwards transform the dimensionless 
free parameters back into the physical variables that char- 
acterize the system. This procedure ensures that its outcome 
contains all the correlations that exist among the physical 
parameters of dusty winds. 

Our presentation starts with a single type of dust grains, 
section |^ extends the discussion to mixtures of sizes and 
chemical compositions. The grain is specified by its size a, 
condensation temperature Tc and absorption and scatter- 
ing efficiencies Qabs,A and Qsca,A. We associate the radius 
ri with prompt dust formation so that Td('"i) ~ Tc- The 
mathematical problem does not contain any reference to ri , 
its actual magnitude enters only during the final transfor- 
mation to physical quantities. We address this issue again in 
our summary in section ^ Table 1 lists the dust properties 
used in our numerical applications. Appendix ^ contains a 
glossary of all the relevant symbols. 



2.2 Dynamics 

We start with a discussion of the three force components 
controlling the supersonic phase, identifying in each case the 
characteristic velocity scale and the dimensionless profile of 
its radial variation. These physical processes have been dis- 
cussed extensively in the literature, most recently by NE 
and HTT. We repeat the essential ingredients to establish 
the proper formalism and examine the various assumptions 
underlying the theory. We also offer a slightly improved ex- 
pression for the dust drift velocity, accounting for the sub- 
sonic regime. 



2.2.1 Radiation Pressure 

The radiation pressure force per unit volume is 
1 . , . . 



Here Qpr is the radiation pressure efficiency (= Qabs + Qsca, 
assuming isotropic scattering) and F\ is the local radiative 
flux, comprised of the attenuated-stellar and diffuse contri- 
butions; note that the diffuse flux vanishes at ri (IE97). The 



spectral matching between F\ and the dust opacity varies 
in the wind because of the reddening of the radiation. This 
variation is conveniently described by the following radial 
profile, normalized to unity a.t y = 1: 



Here F = J F\ d\ is the bolometric flux and 

Qpr,xBx{n)d\, 



}pr,X J^!'.. } dX ■ 



F{1) 



(3) 



(4) 



where B\ is the Planck function. Table 1 lists the values of 
Q* for our standard grains and T* = 2500 K; the dependence 
on T* is insigniflcant when this quantity is varied within its 
physical range. Because of the reddening of the radiation, 
the spectral matching tends to decrease with radial distance 
so that 0(2/) ^ 1. 

The velocity scale associated with the radiation pressure 
force is defined via Vp = 2 ri.7^i.ad('"i)/p(f'i). Introduce 



2 f^d 
2 



= 10 



-22 2 

(722 cm , 



(5) 



the cross-section area per gas particle upon dust condensa- 
tion. Then 



^27rmpCri 



1/2 



111 km s" 



Qi,a22Li 
ri,i4 



1/2 



(6) 



where 7-1,14 = ri/lO^'* cm. Ignoring momentarily the drag 
and gravity effects, the radiative force gives the equation of 
motion 

2 

(7) 



dv' v~ 



When reddening is neglected too, cf> = 1 and the solution is 
simply 



2 2.2(1 

V = «T + «p 1 

' y 



(8) 



Here we take as the starting point for the outflow velocity 
the isothermal sound speed 



TTlH, 



1/2 



= 2.03 T^^^ km s" 



(9) 



where Tjj = T^a x 1000 K is the kinetic temperature at y 
= 1 (not necessarily equal to the dust temperature at that 
point). In this approximation, first derived in GS, the out- 



flow flnal velocity is {vq, + Vp 



2x1/2 



since vt is usually 



negligible. Typical values for the free parameters produce 
a velocity scale Vp considerably higher than observed out- 
flow velocities, a problem noted by Castor (1981) as a se- 
rious shortcoming of dust-driven wind models. The effects 
of reddening, drift and gravity must supplement radiation 
pressure for a viable explanation. 



2.2.2 Drag 

Collisional coupling accelerates the gas particles and decel- 
erates the dust. Gilman (1972) has shown that the dust- 
gas relative velocity reaches steady state within a distance 
£ <^ ri and that the dust then fully mediates to the gas the 
radiation pressure force. In steady-state drift the drag force 
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can be eliminated and the separate equations of motion for 
the dust and gas combined into a single equation. HTT pro- 
vide a useful discussion of this approach, which is the one 
taken here. Still, the dust drift has an important effect be- 
cause the radiative acceleration is proportional to nd/nn, 
and separate mass conservation for the dust and the gas im- 
plies Wd/riH oc v/vd; in spite of our assumption of prompt 
dust formation and no further formation or destruction, the 
dust abundance varies in the shell because of the difference 
between the dust and gas velocities. 

In appendix ^ we derive a simple expression for the drift 
velocity, including both the subsonic and supersonic regimes. 
The drift effect introduces the independent velocity scale 



Mc 



203 km s" 



Q.L4 

M-6 



and the steady-state drift velocity is 



(10) 



(11) 



(equation B4). The significance of subsonic drift is rapidly 



diminished with radial distance because v is increasing as the 
gas accelerates while vt is decreasing as it cools down. The 
radial variation of vt requires the gas temperature profile, a 
quantity that does not impact any other aspect of the flow 
and whose calculation contains large uncertainties. We avoid 
these uncertainties and use instead the initial vt throughout 
the outflow. This slightly overestimates the overall impact 
of subsonic drift, producing a negligible error in an effect 
that is small to begin with. 

The dust velocity is Vd ~ v+Vid, therefore Ud/nn varies 
in proportion to the dimensionless drift profile 



ay) = - 



(12) 



Since na^nd{y)/nH{y) ~ o"gC(y) for J/ li the force equation 
including the drift effect is 



dy y^ 



(13) 



Because of the drift, at y = 1 the radiative force is reduced 
by a factor C(l)i ^ significant reduction when Um S> wt. The 
reason is that the dust particles are produced at the veloc- 
ity VT with a certain abundance and the drift immediately 
dilutes that abundance within a distance ^ ^ ri so that 
nd/nn is diminished already at y = 1. In calculations that 
keep track separately of the dust and the gas, the two species 
start with the same velocity and this dilution is generated 
automatically (cf NE, HTT). Here we solve for only one 
component, the two species start with different velocities at 
y = 1 and the initial dilution must be inserted explicitly. 

It is important to note that C is a monotonically in- 
creasing function of y (see figure y below). Therefore, the 
dust abundance is the smallest at y = 1 and increases from 
this minimum during the outflow. At the wind outer regions, 
rid/nn exceeds its initial value by the factor C(oo)/C(l). 



2.2.3 Gravity 

The gravitational pull introduces an independent velocity 
scale, the escape velocity at ri 



(^) 



15.2 km s" 



Mq 

ri.i4 



1/2 



(14) 



which is frequently entered in terms of the dimensionless 
ratio 



r = 



^ 2:1 



= 45.8(5,(722 



Mo 



(15) 



Gravity does not introduce any radial profile beyond its 
variation. Adding the gravitational effects, the equation of 
motion becomes 



Hy)ay) - p 



(16) 



dv^ Vp 
dy 1/2 

In the limit of negligible drift and reddening ((^ = = 1), the 
gravitational pull reduces the wind terminal velocity from Vp 
to Vp{l - 1/r)^/^ typically only ~ 1% effect. 

This is the complete form of the equation of motion, in- 
cluding all the dynamical processes in the wind. The outflow 
is fully specified by the four independent velocity scales Dp, 
Hm, '^g and Vt, which is also the initial velocity v{y = 1), and 
the reddening profile (p. This profile varies with the overall 
optical depth and is determined from an independent equa- 
tion, the equation of radiative transfer. 

2.3 Radiative Transfer 

Because of the spherical symmetry, the radiative transfer 
equation requires as input only two additional quantities 
(IE97). One is the overall optical depth along a radial ray 
at one wavelength, which we take as visual 



TV = ricrv / riddy 



(17) 



the optical depth at every other wavelength is simply 
TvQ\/Qv, where Qx and Qv are the absorption efficiencies 
at A and visual, respectively. The other required input is the 
normalized radial profile of the dust density distribution 

ridiy) 



viy) = 



(18) 



Uddy 



Given these two quantities, the intensity Ix(y, (3) of radiation 
propagating at angle /3 to the radius vector at distance y can 
be obtained from 

dIx{y,P) _ ^ , 



drx{y,f3) 
where 

T\{y,l3) = TV 



^ cos /3 



(19) 



(20) 



In general, 77 and ry must be specified as independent input. 
Instead, here they are fully specified by the dynamics prob- 
lem. With the assumption of no grain growth or sputtering 
after its prompt formation, mass conservation for the dust 
gives Ud oc 1/r^Vd oc C,/y'^v so that 



■n{y) 



gy) 

y^v{y) 



Ldy 



y-v 



(21) 



Therefore -q is not an independent input property, instead 
it is determined by the solution of the equation of motion. 
And straightforward manipulations show that 
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J, y-'v 



(22) 



so TV, too, does not require the introduction of any addi- 
tional independent input. The radiative transfer problem 
is fully specified by the parameters already introduced to 
formulate the equation of motion; it requires no additional 
input. 

This completes the formulation of the dusty wind prob- 
lem. The dynamics and radiative transfer problems are cou- 
pled through the reddening profile (f). The solution of the 
equation of motion obtained with a certain such profile de- 
termines a dust distribution and an optical depth rv- The 
solution of the radiative transfer equation obtained with 
these rj and ry as input properties must reproduce the red- 
dening profile that was used in the equation of motion to 
derive them. 



2.4 Scaling 

Given grain properties and the spectral profile of the stel- 
lar radiation, the wind problem is fully specified by the four 
independent velocities Vp, Vm, Vg and vt, which together 
form the complete input of the problem. Since an arbitrary 
velocity magnitude can always be scaled out, the mathemat- 
ical model can be described in a dimensionless form which 
includes only three parameters. We choose Vm for this pur- 
pose because it results in a particularly simple equation of 
motion. Introduce 



w = v/i 



' = vt/v„ 



P = Vp/Vn 



(23) 



The parameter P characterizes the ratio of radiation pres- 
sure to drift effects. It is very large when the drift becomes 
negligible and goes to zero when the drift dominates. The 
equation of motion ( |l^ ) becomes 



dy 



(0C 



where C, ■ 



9 + ^/w4> 



(24) 



and where is determined from the solution of the radiative 
transfer equation (^^ in which the dust distribution and 
optical depth are 



v{y) = 



y'^w{y) 



TV 



QvP' 



^— dy 

y^w 

y^w 



(25) 



In this new form, the mathematical model is fully specified 
by the three dimensionless parameters P, T and 9, where 
the latter is also the initial value 9 = w{y = 1). The velocity 
must rise at the origin, and since — 1 the parameters 
are subject to the constraint 



rc(i) > 1, 



(26) 



This liftoff condition ensures that radiation pressure, with 
proper accounting for the dust drift, can overcome the initial 
gravitational pull. It can be viewed as either a lower limit 
on r for a given 6 or, given F, a lower limit on 6. 

The dusty wind problem has been transformed into a 
set of two general mathematical equations — the equation 
of motion ( p^ and the radiative transfer equation ( |l9[ ) with 
the input properties from equation (ESI). Given the spectral 



shapes of the stellar radiation and the dust absorption co- 
efficient, these equations are fully prescribed by P, T and 
9. Both the radiative transfer and the dynamics equations 
are solved without any reference to Vm or any other velocity 
scale. The velocity scale is an extraneous parameter, entirely 
arbitrary as far as the solution is concerned. The wind ra- 
diative emission and the shape of its velocity profile are both 
independent of the actual magnitude of the velocity. In the 
following discussion we show that the shape of w{y) turns 
out to be a nearly universal function, independent of input 
parameters within their relevant range, and that the final 
velocity Wac = w{y ^ oo) is for all practical purposes only 
a function of P. 



3 SOLUTIONS 

The complete solution involves two elements — dynamics 
and radiative transfer. The impact of the radiative transfer 
on the dynamics can be conveniently expressed in terms of 
the quantity $, defined via 



Woe 



dw 



(27) 



That is, $ is the velocity-weighted harmonic average of the 
reddening profile 0. Combining equations ^ and^^, it is 
easy to show that 



+ 



2r 



dy 

y^w(j) 



(28) 



Therefore, when gravity is negligible (F ^ 1) the optical 
depth and the terminal velocity obey the simple relation 



Qv 



(29) 



In optically thin winds <I> = 1 and Woa = tv/Qv- As red- 
dening increases $ decreases, the final velocity is affected 
by the radiative transfer and numerical computations are 
necessary. 

We obtain numerical solutions for arbitrary optical 
depths with the code DUSTY (Ivezic, Nenkova & Elitzur, 
1999). DUSTY solves the dusty wind problem through a 
full treatment of the radiation field, including scattering, 
absorption and emission by the dust, coupled to the hydro- 
dynamics problem formulated above. Appendix ^ provides 
a description of the numerical procedure. We present now 
the solutions, starting with the optically thin regime where 
we found the exact analytic solution. 

3.1 Negligible Reddening 

When reddening is negligible, the flux spectral shape does 
not vary and </!> = 1 (see eq. This is the situation in 
optically thin winds, and the detailed results presented in 
the next section show this to be the case when rv ^ 1.Q 
Even though the optical depth at wavelengths shorter than 
visual already exceeds unity when rv = 1 and the emerging 

* Another limit that in principle could give </> = 1 involves grains 
so large that Qx is constant for all relevant wavelengths. This 
requires grain sizes in excess of 10 |am, and seems of little 
relevance. 
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spectrum is affected by radiative transfer, the impact on 
the dynamics is minimal. The reason is that (j> involves a 
spectral average, and the stellar radiation peaks at longer 
wavelengths. 

The wind problem with tf) — 1 was called the "first sim- 
plified model" by HTT, who solved it numerically. Here we 
present the complete analytic solution for this case. Since (j) 
is known, the equation of motion decouples from the radia- 
tive transfer problem and can be considered independently. 
Introduce 



(30) 



r- 1' 

then equation ^ can be cast in the form 

1 + 5 \ dw^ _ 1 
w + e-5 ) ly ^ 1 + 5 ^ 



1 + 



so that 



- e"^ + A{1 + 5) 



X dx 

x + e- 



1 + 5 



1 - 



(31) 



(32) 



The liftoff constraint ( Pq) ensures that the denominator of 
the integrand never vanishes in physical solutions; it starts 
positive at the lower boundary and increases as x increases. 
The integration is standard, the result is 



1 + 5 



1 - 



= w^- 61^+4(1-^5) 



+ \{5-e)(w-> 

1/2 



i(«^/^ 



+ (5 -eY(w''^ -e^''^) 
+ [5-ef\n'^^^-^ 



n3/2 



Ve + e-s 



(33) 



This is the complete solution for all dusty winds with rv < 1, 
fully incorporating the effects of gravity and dust drift. Radi- 
ation reddening, which takes effect when ry > 1, is the only 
ingredient missing from this analytic solution and preventing 
it from applicability for all winds. The GS solution (equa- 
tion ^), the only previous analytic result, which neglected 
the effects of gravity and drift in addition to reddening, can 
be recovered inserting 5 = and P ^ 1. 

The solutions are meaningful only when describing 
winds in which Voo/vt is at least ~ 3; otherwise, the effects 
of gas pressure, which were neglected here, become signifi- 
cant and equation ^ loses its physical relevance. Therefore, 
we require «oo/«t ~ Woo/d > 3. In addition, the liftoff con- 
dition implies that 5 < 5max where 

<5max = Ve + e. (34) 

Figure |l| displays the dimensionless velocity profile as a func- 
tion of distance y for a range of representative values of 
the three free parameters and demonstrates a remarkable 
property: Except for their role in determining the boundary 
of allowed phase space and controlling the wind properties 
near that boundary, 5 and 6 hardly matter; away from the 
boundary, the solution is controlled almost exclusively by 
the single parameter P. Part (a) of the figure shows the ef- 
fect of varying 5 when P and 6 are held fixed. Each panel 
presents a representative value of P with a wind initial ve- 
locity 6 — 0.03P; this choice of 9 ensures that the wind 



velocity increase by a factor of ~ 10-20 in each cases. The 
displayed values of 5 cover the entire physical range for this 
parameter, from to just below the singularity at Jmax- In 
each panel, the solutions for 5/(5max ^5 5 are rather similar 
to each other, those with 5/5max ^5 j hardly distinguishable. 
All display a similar rapid rise within y ^ 10 toward a fi- 
nal, nearly the same velocity. The plots for 5/(5niax 0.999 
stand out and show t he q uenching effect of gravity, which is 
discussed below (sec. 5.1). 

The effect of 6 on the solution when 5 and P are held 
fixed are similarly shown in part (b) of the figure. In these 
panels 5 = 0. This value was chosen because it represents 
faithfully most non-quenched solutions while allowing the 
wind to start with arbitrarily small initial velocity, even zero. 
The values of 9 for the displayed solutions are listed as frac- 
tions of Smax, the initial velocity that leads to Woa/d = 3 
for the listed P. As is evident from the plots, the initial ve- 
locity is largely irrelevant. Starting from S = or as much 
as |6'max yields practically the same results. The only dis- 
cernible difference occurs at ^max, where the profile deviates 
from the common shape by no more than ~ 10%. 

The analytic solution defines w as an implicit function 
of y. In appendix |c| we derive explicit analytic expressions 
for w that provide adequate approximations in all regions of 
interest. An inspection of the solution shows that P = 16/9 
is the transition between the drift-dominated regime at small 
P and negligible drift at large P. The velocity profiles in the 
two are 



W — Woo I 1 

y 



2/3 for P < 16/9 
1/2 for P> 16/9 



(35) 



The expressions for Woo in the two regimes can be combined 
into the single form 



= P 



P 



1/3 



(36) 



The optical depth, needed for the radiative transfer prob- 
lem which is solved independently in this case, is simply rv 
= QyWoo (see eq. p9[ ) . Appendix |c| provides the 5- and 9- 
corrections to these results, which show that the dependence 
on 9 is confined to the very origin of the wind, y — 1^ 6/woo ■ 
The last two expressions reproduce to within 30% the 
plots for all the non-quenched winds displayed in figure |l[ 
They amplify our conclusion about the negligible role of 
both 5 and 6 away from the physical boundary and point 
to another remarkable property of the solution. Although 
the final velocity strongly depends on P, the shape of the 
velocity profile w/woo does not. The only reference to P in 
this profile is in determining a transition to a slightly steeper 
profile shape at large P. The profile itself is independent of 
P in either regime and furthermore, the difference between 
the two shapes is not that large. These properties are evi- 
dent from figure ^. Other than the numerical values on the 
vertical axes, it is hard to tell apart the plots in the different 
panels. 

These results show that among the three independent 
input parameters, P is the only one to have a significant ef- 
fect on any property of interest. Furthermore, even P affects 
mostly just the final velocity, it does not have a discernible 
effect on the shape of the velocity profile; all optically thin 
winds share universal velocity and dust density profiles. 
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(a) 



(b) 



0.01 -i 




0.01 



1000 




Figure 1. The complete solution (equation p3| ) for the velocity profiles of optically thin winds, (a) In each panel P has the indicated value 
and the initial velocity is 9 = 0.03P. The solutions are plotted for various values of 5, marked as fractions of the gravitational-quenching 
value (5max. Note the virtual 5-independence, except for the solutions asymptotically close to 5max. (b) Same as (a), only now 8 varies 
and S is fixed at 0. The values of 6 are listed as fractions of 6max, the initial velocity that yields the smallest meaningful velocity increase 
Wao/0 = Voo/vt = 3. Note the complete 9-independence, except for the solution with 9 = 0max- 



3.2 Reddening Effects 

The results for optically thin winds carry over to arbitrary 
optical depths, with profound implications for all winds. 
Reddening effects are controlled by the wind optical depth 
and density profile. For a given pair of S and 0, consider a 
value of P sufficiently small that rv ^ 1 so that redden- 
ing can be neglected; such a choice of P is virtually always 
possible. The velocity profile and rv are then uniquely de- 
termined by P (equations ^ and [s^ ) and so is the radiative 
transfer problem. When P increases, ry increases too. The 
results of Appendix O show that 




yields ry ~ 1; with our standard grains, the corresponding 
values are P — 0.73 for carbon and P = 1.35 for silicates. As 
P increases further the optical depth increases too and with 
it the impact of reddening, and that impact too is controlled 



exclusively by P. Therefore, the parameter P controls all 
aspects of the problem. As concluded in 1E95, away from 
its boundary and for most of phase space the dusty mnd 
problem is controlled by a single free parameter. Since P, tv 
and Woo are uniquely related to each other, either one can 
serve as that free parameter. 

Figure |^ shows the radial variation of profiles of inter- 
est for various values of P. The models for P — 0.1, 1 and 
10 repeat those presented in figure |l] but fully incorporate 
radiative transfer. Comparison of the corresponding plots in 
the two figures illustrates the impact of reddening. The bot- 
tom panel of figure ^ shows the reddening profile cj)- Most 
of the reddening occurs close to wind origin. Reddening be- 
comes significant at P = 1, and substantial at P = 10 and 
100. Still it has only a minimal impact on the velocity pro- 
file; as is evident from the top panel, the dependence of the 
profile shape on P remains weak. We find that equation ^ 
remains an excellent approximation under all circumstances. 
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Figure 2. Radial profiles for quantities of interest at various val- 
ues of P, as marked. The corresponding optical depths are shown 
in parentheses. All models have amorphous carbon grains, 5 = 
and Woc/6 = 10. 



the effect of reddening merely decreases the exponent k from 
0.5 to 0.4 at large P. The plots of show that most of the 
drift variation, too, occurs close to the origin. At y ^ 2 the 
dust and gas velocities maintain a constant ratio, which can 
be substantial when P is small; at P = 0.1 the final veloc- 
ity of the dust grains exceeds that of the gas particles by 
more than factor 6. The velocity and drift profiles combine 
to determine the dust density profile rj (see eq. ESI). The fig- 



Figure 3. Variation with P of various quantities of interest. 



ure presents the fijnction y r], removing the common radial 
divergence factor 

The plots present the full solutions of the dynamics 
problem coupled with radiative transfer for carbon grains. 
The results for silicate grains are essentially the same, ex- 
tending the conclusion reached in the optically thin limit 
to all cases: all radiatively driven dusty winds share nearly 
identical velocity and dust density profiles. 

Figure ^presents the dependence of tv/Qv {~ Woo/^) 
and final velocity vUoo on P. The bottom panel shows the 
reddening indicator <1>. As long as <1> = 1, reddening is neg- 
ligible, the analytic solution (B3) holds and the grain prop- 
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erties are irrelevant. Indeed, all quantities are the same for 
carbon and silicate dust in that regime. The grain composi- 
tion matters only when reddening affects the dynamics. The 
figure shows that ^ begins to decrease, indicating that red- 
dening becomes significant, at a value of P corresponding to 
TV — I; the actual value is controlled by Qv (equation ply , 
it is larger for silicate because of its smaller Qv. As noted 
before, the optical depth at wavelengths shorter than visual 
already exceeds unity when rv = 1, but that spectral region 
has little impact on the dynamics. 

Appendix B shows that the analytic approximation in 
equation ^ can be extended to include reddening effects, 

/ P V'" 
Woo = pV^ TT-^ ■ (38) 

This expression provides an excellent approximation for the 
numerical results presented in figure pi the deviations are 
less than 10% for silicate grains, 30% for carbon. In addition 
to Waa, the figure displays also Wao/P = Voa/vp. As discussed 
in ^ the simplest application of radiation pressure gives the 
GS result Voo = Vp, i.e. Woo/P ~ 1- The figure shows that 
the actual solution has a substantially different behavior — 
Woo/ P is neither constant, nor does it reach unity. Instead, 
it increases in proportion to P^/^ at P ^ 1, because of the 
drift, and decreases in proportion to ^^^^ at P ^ 1, because 
of the reddening. Its maximum, reached around P ~ 1, is 
larger for silicate because it requires larger P for reddening 
to become significant. Drift and reddening play a crucial 
role in the behavior of Woo/P- The maximum reached by 
this function translates to a maximum velocity for dusty 
winds, discussed in the next section. 



4 OBSERVABLE CORRELATIONS 

Contact with observations is made by transforming the 
mathematical problem back into physical quantities. Except 
for regions close to the boundary of phase space (see the 
wind problem is fully controlled by the parameter P and its 
solution determines the corresponding Woc- Both quantities 
will now be expressed in terms of physical parameters. 
The main scaling variable is 

1/2 

P = 0.546 A/_6 ( ^ I (39) 



Q, 1/4 ri,i4 

(see equation |23[ ) . This expression does not determine P ex- 
plicitly because P enters indirectly also on its right hand 
side. The dust condensation condition T{ri) = Tc deter- 
mines ri as (IE97) 



ri = 



1/2 



1.16x10' 



rp2 

^c3 



(40) 



Here Tc3 — Tc/(1000 K) and is a dimensionless function 
determined by the radiative transfer, similar to the redden- 
ing profile, and thus dependent on P. Optically thin winds 
have 'I' = "Jo where 

Here Qp{T) is the Planck average of the absorption effi- 
ciency, similar to the average of the radiation pressure ef- 
ficiency that defines Q, (equation W); Table 1 lists 'I'o for 



our standard grains. As P increases, the wind becomes op- 
tically thick and 5' increases too. In IE97 we present the 
variation of ^/'I'o with optical depth and show that it is 
well approximated by the analytic expression 



* =i + 3^Qp(r.; 

where 



(42) 



(see figure 1 and equation B7 in IE97)Q Just as <I> is the 
velocity-weighted harmonic average of the reddening profile, 
(f> is its standard average, weighted by r]/y^. 

For most of the relevant region of phase space, P fully 
controls the solution of the dusty wind problem. Since ^' is 
part of that solution, it too depends only on P. Therefore, 
the combination of equations M and klOl results in 



,1/2 



r3/4 



= 1.! 



(Too Jc3 



P* 



1/4 



(43) 



a one-to-one correspondence between M/L^^'^ and P. This 
correspondence involves the function P^^^'*, which is fully 
determined from the solution of the mathematical wind 
problem. The proportionality constant involves the individ- 
ual grain properties Q, and Tc which are part of the prob- 
lem specification, and (J22 which is not. The dust abundance, 
necessary for specifying CT22, has no effect on the wind solu- 
tion; only the proportionality constant is modified when this 
abundance is varied, the function Pvl'^'''' remains the same. 

This result shows that for any dusty wind, the combi- 
nation M/L'^^* can be determined directly from the shape 
of the spectral energy distribution: the wind IR signature 
is fully controlled by the parameter P, therefore comparing 
the observations with a bank of solutions determines the 
value of P and equation ^ fixes M/L^^^. We may expect 
all C-rich stars to have dust with roughly similar abundance 
and individual grain properties, so that they have the same 
proportionality constant and functional dependence ^'(P); 
likewise for O-rich stars. In that case, each family defines a 
unique correspondence between M/L^^* and P. And since 
the parameter P fully controls the IR signature, stars that 
have different M and L but the same Af/L"^^* are indistin- 
guishable by thetr IR emission because they also have the 
same P. 

The solution of the wind problem also determines 
Woo, another unique function of P. Therefore, the relation 
Voo/vm = Woo gives another one-to-one correspondence be- 
tween P and an independent combination of physical quan- 
tities 



L 



(44) 



Since this correspondence does not involve o"22 , it is indepen- 
dent of the dust abundance and depends only on individual 
grain properties. Systems with different M , L and Voo but 

the same M jL^^'^ necessarily have also the same Mvoo/L 



t Since scattering is neglected in the IE97 analytic solution, the 
difference between Q« and Qp(T«) is ignored in the approxima- 
tions here (but not in the numerical calculations). 



© 2001 RAS, MNRAS 000, [|-|^ 



10 Moshe Elitzur & Zeljko Ivezic 



because each combination uniquely determines the parame- 
ter P of the system. 

The complete description of the wind problem involves 
four velocity scales, two of which {vt and Vg) do not count 
whenever F and 9 can be ignored. For most of phase space 
the problem is fully specified by the velocities Vp and Vm, 
and its solution determines Voa ■ The two ratios formed out of 
these three velocities are all the independent dimensionless 
combinations of physical parameters to determine the math- 
ematical wind model and its solution. Since this solution 
fully specifies the wind IR signature, M/L^^* and Mvoo/L 
are the only combinations of global parameters that can be 
determined from IR observations, even the most detailed 
ones. When the velocity is additionally measured in molecu- 
laJT line observations so that both P and Voo are known, the 
two combinations in equations ^ and ^ can be used to de- 
termine also L and M individually. The relevant correlations 



1^1/4 p^l/4 



(45) 



whose constants of proportionality can be trivially derived. 
Only two of the last four combinations that relate M, L and 
Voo to the solution are independent; any two of them can be 
derived from the other two. 



4.1 Similarity Relations 

While P is the natural independent variable of the mathe- 
matical problem, its physical interpretation is only indirect. 
Modeling of IR observations directly determines rv, not P. 
Therefore, we now express all results in terms of ry by per- 
forming a straightforward change of variables. The relation 
between rv and P is non-linear, starting from P oc Ty* in 
the optically-thin regime and switching to a more complex 
dependence when reddening becomes significant. We intro- 
duce the Ty-P transformation function Q through 



^/3 Wv ' 



(46) 



so that = 1 when rv < 1. At larger optical depths, B 
can only be determined from the numerical solution. The 
analytic expression 



[l + Krv/Qv)'/']*, 



(47) 



obtained from equation C12, provides a useful approxima- 
tion for the numerical results at all optical depths. With the 
Ty-P transformation, equation ^ becomes 



M- 



j-3/4 

where 



= Ci Ty^^Ki{Ty), 



(48) 



l/2,T,l/4 



ci = 2.28 



^3/4 l/2rp 

^y ^22 -'cS 



From its definition K\ — 1 when rv < 1, therefore M_6 = 
ci(L4rv)'^''* in all optically thin winds. As the optical depth 
increases, Ki introduces the reddening correction to this 
relation. This correction, purely a function of rv, is shown 
in the top panel of figure H. 



The transformation of the other independent correla- 
tion from P to Ty is considerably simpler. Thanks to equa- 
tion ^ which holds for all winds when gravity is negligible, 
equation ^ can be cast similarly in the form 



Mvo 



■iQ,/Qy)TyK2iry), 



(49) 



where K2 = $. The function K2, shown in the second panel 
of figure 1^, contains the reddening corrections to the opti- 
cally thin correlation, similar to the previous result. Since 
K2 = $, it is the same function as shown in figure |^, only 
plotted against the independent variable rv instead of P. 

This completes the transformation from P to ry of 
the two independent correlations. In both, the luminos- 
ity does not enter independently, only through the prod- 
uct Lry. In the case of equation ^ this is a direct conse- 
quence of the structure of the equation of motion, which 
gives = v^Wao = r}^Ty{<^/Qy). The combination Lry 
emerges here as a "kinematic" result irrespective of the ex- 
plicit expression for the drift function. In the case of equation 
p^ , on the other hand, this combination is a direct conse- 
quence of the specific functional form of C,. With these two 
results, the correlations in eq. ^ become similarly 



^1 

7-1/4 



■ C3 ry''_K"3(rv), 



C4 A'4(rv) 



(50) 



where K-^ = K2/K1 and — K2/KI. The associated pro- 
portionality constants are 

C3 = 88.9rc3(O.O-22)''''(Qv*0)"'''^ 

C4 = 3.08X 10^ T.ljQ.aia*-! . (51) 

Since L enters only in the form Lry, the combination 
that eliminates L eliminates also ry from its optically thin 
regime. The reddening correction functions and K4, are 
shown in the two bottom panels of figure ^ 

At small Ty the grain properties are irrelevant and the 
similarity Jf-functions are unity for both silicate and carbon 
dust. As Ty increases the curves for the two species diverge, 
refiecting their different optical properties. However, figure 
^ shows that the differences are quite moderate in all four 
cases. They are most noticeable in K\ and K2, where the 
relative differences from the mean never exceed ~ 25%. It 
is also important to note that the displayed range of rv 
greatly exceeds the observed range since cases of rv <; 10 
are rare. The differences are greatly reduced in the ratio 
function (= K2/K1) and practically disappear in K4. 
This function is essentially the same (within a few percent) 
for both silicate and carbon grains at all optical depths, 
an agreement maintained over many orders of magnitude. 
Evidently, K4 is independent of chemical composition. 

The universality of K4 is no accident. From the ana- 
lytic approximations for (eq. ^ ) and O (eq. ^ ), at large 
optical depths TyKi cc ^/4>\ that is, apart from the explicit 
dependence on rv, the variation of Ki comes from the ratio 
of two averages of the reddening profile, one ($) strongly 
weighted toward the outer part of the shell the other (<^) 
toward the inner part. Although $ and 4> are different for 
grains with different optical properties, their ratio is con- 
trolled by the shapes of the density and velocity profiles — 
which are universal, as shown above. As a result, the func- 
tional dependence of K4, on rv is the same for all grains 
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Figure 4. The reddening corrections in the scaUng relations 
that summarize the physical contents of the solution (equations 
^8|-po|). Note the linear scale in the plot of Ki. The full and 
dashed lines show the results of detailed numerical calculations 
for carbon and silicate grains, respectively. The dottcd-line in the 
panel for Ki (« = 1 . . . 4) plots the function 1/(1 + tv)°' where 
ai = -0.105, 02 = 0.36, as = 0.465 and C14 = 1.5. 



when TV S> 1 in addition to tv < 1. This does not yet guar- 
anty a universal profile because the transition between the 
two regimes could depend on the grain parameters. Indeed, 
the plot of $ {— K2) shows that the large- rv decline of this 
function is delayed for silicate in comparison with carbon 
dust. As explained in the discussion of figure 0, the onset 



of reddening requires a larger rv for silicate because of its 
smaller Qv. But the behavior of the factor 1 -j- |('rv/Qv)^''^ 
in O (eq. is precisely the opposite since its large- rv be- 
havior starts earlier when Qv is smaller. The two effects 
offset each other, producing a universal shape. 

The simple analytic expression {l + T-\f)~^'^ fits the func- 
tion K4, to within 20% over a variation range covering more 
than four orders of magnitude, providing the nearly perfect 
fit evident in figure ^. Therefore, dusty winds obey the con- 
dition 



M- 



Vf , ^1.5 

(1 +rv) 



C4 . 



(52) 



This result provides a complete separation of the system 
global parameters from the grain parameters. The combi- 
nation of Voo, M and rv on the left hand side is always 
constant, its magnitude is determined purely by the dust 
properties. This expression makes it evident again that the 
optically thin regime corresponds to rv < 1. 

As is evident from figure K3 too is nearly the same 
for carbon and silicate grains. The reason is that K3 = 
{K2K4y^'^ and the only difference between the two species 
comes from their K2 profiles, which enter only in the fourth 
root. The single analytic expression (1 -j- rv)"" *®^ provides 
an excellent fit for both silicate and carbon grains, leading 
to the independent correlation 



ui = C3 (L4rv)^'''* (1 + rv)" 



(53) 



This result can explain the narrow range of velocities ob- 
served in dusty winds. The dependence of velocity on rv 
at fixed luminosity and grain properties is shown in figure 
^ This function reaches a maximum of 0.73 at rv = 1.3, 
therefore the largest velocity a dusty wind can have is 



= 0.73 C3 L^^ 



km I 



(54) 



Figure |^ shows that the deviations from this maximum are 
no more than a factor of ~ 2 when rv is varied in either 
direction by two orders of magnitude. The dependence of 
velocity on luminosity is only L^^'^, and since L4 is typically 
~ 0.3-20 it introduces a similarly small variation. Finally, 
even the dependence on grain properties is weak, as is ev- 
ident from the expression for C3 (eq. [si] ). The only grain 
parameter that enters linearly is the condensation tempera- 
ture, and its dispersion is expected to be small. 

Since only two of the _ftr-functions are independent, the 
universal fits for and K4 can be used to produce grain- 
independent approximations also for Ki and K2- Figure ^ 
shows that the single function K^/Ka, = (!-!- tyY''^^^ pro- 
vides an adequate fit for Ki of both silicate and carbon dust, 
always within 30% of the detailed result^. This expression 
gives 



= ci (L4rv)^'"' (1 + rv 



(55) 



Since rv = 1.3 gives the largest possible velocity at a given 
L, the corresponding mass loss rate is 



M-6(lJmax) 



1.33 ci L\'^. 



(56) 



t Higher accuracy, when desired, can be obtained by replacing 
the common index 0.105 with 0.05 for amorphous carbon and 
0.15 for silicate grains. 
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Figure 5. The ry-variation of the correlations jHs] and 



A common practice in analysis of observations is to deduce 
the optical depth from spectral data and derive M from rv 
(e.g., Jura 1991). This procedure is predicated on the as- 
sumption of a linear relationship between M and ry, and 
our result shows that this assumption is not quite right. At 
a fixed luminosity, M is proportional to t^^* in the optically 
thin regime and to Ty^^ at large ry. The expectation of a 
linear relationship between M and tv is not met because 
of the dust drift, since ry/Moc J {(,/y'^v)dy. The relation- 
ship would be linear if were 1 everywhere, but figure ^ 
shows that this is never the case. The deviations from unity 
are large at small ry, where the drift is most prominent, 
decreasing as ry increases. 

Finally, figure W shows that the universal profile 
K^/Ki — (l+ry)"" again describes reasonably the actual 
K2 of both grains, so that 

L 



■{Q,/Qv) rv(l + ry)- 



(57) 



The relation Mv^o ^ L/c has often been used as a physical 
bound on radiatively driven winds, even though the mistake 
in this application when ry > 1 has been pointed out repeat- 
edly. In IE95 we show that the proper form of momentum 
conservation is Mvao = rpL/c where rp is the fiux-averaged 
optical depth, so now we have found the explicit expression 



TP = (Q,/Qy) Tv{l+Tv) 



(58) 



The Ty-variation of this function is shown in figure y. It 
increases linearly at small ry, reaches unity at rv ~ 1.6 
and switches to Ty^* thereafter. The ratio Q,/Qv is 0.1 for 
silicate and 0.25 for carbon (see Table 1). 

For every pair among M, L and Voa, equations (js^), ( |53| ) 
and ( |55| ) list the correlation in term of optical depth. When 
detailed IR data are not available and ry cannot be deter- 
mined, it can be bypassed by correlating the pairs directly 
against each other. Since only two of the three relations are 
independent, there is only one such combination. Equation 
is the most suitable for eliminating rv because Ki varies 
the least among the four A'-functions. In fact, the crude ap- 



proximation A'l ~ 1 introduces an error of less then 50% for 
ry ^ 10. With this approximation, ry oc M^'^^/L. Inserting 
this result in equation |5^ yields 

-1.5 



3 



AM- 



B 



where A = ca and B 



-4/3 



(59) 



This universal correlation 



summarizes our solution for all dusty winds away from the 
boundaries of phase space. The error in this result, intro- 
duced by the approximation Ki = 1, is less than 50% when 
TV ^ 10. When the observational accuracy warrants higher 
precision, corrections can be readily derived. 

This completes the similarity solution of the dusty wind 
problem. Our results amplify the conclusion of IE95, taking 
it a step further: The solution is fully characterized by opti- 
cal depth. The relations among global parameters M, L and 
Uoo involve universal similarity functions of ry, independent 
of chemical composition. The grain properties enter only in 
the proportionality constants of the similarity relations. We 
derived the similarity functions from solutions for carbon 
and silicate grains, whose absorption efficiencies are widely 
different. Since dust spectral features have a negligible ef- 
fect on overall reddening corrections, these functions should 
describe all interstellar grains with reasonable properties. It 
is gratifying that in spite of its great complexity, the dusty 
wind problem can afford such a simple, explicit solution. 

4.2 Young's Correlation 

Young (1995) conducted a survey of nearby Mira variables 
with low mass-loss rates. He finds a clear, strong correlation 
between outflow velocity and mass-loss rate, but indepen- 
dent of luminosity. The correlation can be parametrized as 
M (X v"^, with a = 3.35. Subsequent observations by Knapp 
et al (1998) corroborate Young's results and find a = 2, al- 
though the scatter in their data is consistent with values as 
large as 3. Remarkably, even though the wind is driven by 
radiation pressure, its velocity is independent of luminosity. 

From equation (or |5^) , at small optical depths our so- 
lution gives M oc ii^independent of luminosity, explaining 
the observational findings. The implied a = 3 is consistent 
with all observational results within their errors. This re- 
sult reflects the central role of drift at small mass loss rates. 
When the drift dominates, ( ~ w^^^ oc (Mn/L)^/^. Neglect- 
ing gravity and reverting to physical variables, the equation 
of motion AIm becomes 



(60) 



dr 

where the proportionality constant contains the grain prop- 
erties. The solution gives v^'^ oc {MLY^'^ /r\. And since r\ 
is proportional to L^^^ (eq. ^), the dependence on luminos- 
ity cancels out. If not for the particular L-dependencies of 
ri and the drift, this would not have happened. 

The observed correlation, in particular its luminosity in- 
dependence, directly reflects the specific dependence of the 
drag force on M, L and v. Other forms for the force would 
produce different correlations. For example, when the drift 
is neglected, the combination {MLv)^^^ is replaced by L, 
leading to the GS result Voo = Vp cx {L/riy^^ (eq. |^. To- 
gether with ri oc L^/^ this yields wj^ oc L, as first noted by 
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Jura (1984). This prediction of a correlation between veloc- 
ity and luminosity independent of mass-loss rate is in strong 
conflict with the observations. One could formally eliminate 
the luminosity with the aid of the momentum flux conserva- 
tion Mvoo = rpL/c (see eq. ^ and subsequent discussion) 
to re- write this result as M cc rpu^. While this bears super- 
ficial resemblance to Young's correlation, the observational 
result emerges only if the variation of rp with L and M is 
ignored. By contrast, at rv < 1 equations and M together 
give rp oc TV oc M / L. 

Young's correlation is a direct reflection of the basic 
physics ingredients that went into the model. It demon- 
strates the importance of drift in dusty winds and provides 
strong support for its underlying theory. 



5 PHYSICAL DOMAIN 

A global property not considered thus far is the stellar mass. 
This quantity does not enter into the deflnition of the pa- 
rameter P and thus cannot be determined in general. The 
mass only affects F (see equation |l^, a quantity that can 
vary by orders of magnitude without a discernible effect on 
any observed property. However, F does play an important 
role in determining the parameter range corresponding to 
actual winds. 

The equation of motion ( p^ ) has many mathematical 
solutions but not all of them are physically relevant. By 
example, consider the analytic solution for optically thin 
winds (equation whose logarithmic term becomes sin- 
gular when 5 — (5nmx (see equation The singularity is 
avoided whenever the numerator and denominator are fi- 
nite and have the same sign. Both cases are mathematically 
acceptable but the negative sign is physically meaningless. 
Such solutions violate the liftoff condition (equation ^) and 
the only transition from one set of solutions to the other is 
through the singularity. 

At the outset, P, F and 6 must be positive. Further, all 
winds must obey the liftoff condition, which involves only 
two out of the three input parameters. This condition en- 
sures liftoff under all circumstances but it does not automat- 
ically guaranty a meaningful outcome. The formal solution 
of equation p3 gives 



10 



y 



(61) 



Obviously, in physical solutions the final velocity must ex- 
ceed the initial velocity and we further require Voo/vt 3 
(see §B.l). Because 0^1 and C < 1, {4>C/y'^)dy < 1 and 



P'> 



(62) 



The three parameters of every dusty wind whose velocity 
increases by the ratio Voo/vt must also obey this constraint. 



5.1 Gravitational Quenching 

The wind acceleration is positive as long as F<^(^ > 1 (eq. 
p^l ). Negative acceleration, leading to wind quenching, can 
be caused by increasing either the gravitational pull (smaller 
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Figure 6. Gravitational quenching because of reddening. Models 
are for silicate dust with the listed values of P (the corresponding 
TV is in parenthesis). In each case 9 and F are adjusted so that 
Wca =10 6 and the gravitational pull reaches 99% of its possible 
maximum; the detailed procedure is described in Appendix 



F) or the reddening (smaller (f>). The quenching process de- 
pends on all three input parameters and takes different forms 
in the optically thick and thin regimes. In optically thin 
winds (f> = 1 and since increases monotonically with w (cf 
figure the liftoff condition ensures positive acceleration 
everywhere. Consider the analytic solution (eq. p3[ ) when 5 
increases while P and 6 remain fixed. When /3 = S/Sma.^ ap- 
proaches unity the acceleration becomes negligible, though 
it remains positive, producing the behavior seen in the plots 
for = 0.999 and 0.99999 in figure §. This quenching effect is 
caused by the logarithmic term. Near y = 1 this term domi- 
nates and the velocity increases very slowly. At a distance y 
roughly proportional to P~^6^^^\ ln(l — /3)|, the terms inde- 
pendent of 5 take over and the acceleration picks up. From 
that point on the solution resembles non-quenched winds, 
albeit at a lower acceleration because of its late start. As (5 
gets closer to unity, the dominance of the logarithmic term is 
slowly extended further out until the whole wind is stalled. 

Reddening introduces an entirely different quenching 
mode, affecting winds that obey the liftoff condition when 
P is increased beyond a certain point. The liftoff condition 
ensures that the initial acceleration is always positive, ir- 
respective of optical depth. However, subsequent reddening 
can reduce cj) substantially (figure |^ and since the accelera- 
tion is proportional to cf>C, — l/F, the wind may be prevented 
from reaching a significant terminal velocity. Furthermore, 
since < (j) the acceleration becomes negative whenever 
<j) < 1/T and the velocity can even decline after its initial 
rise. Figure ^ shows examples of such winds whose velocity 
profile differs greatly from the monotonic rise that typifies 
all solutions not too close to the boundary of phase space. 
In each case Woa = 10 8. But when P ^ 25, the velocity 
reaches intermediate values substantially higher than Woa, 
almost as much as ten times higher, before declining to its 
final value. This behavior has intriguing observational impli- 
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cations, but it is not clear whether the limited phase space 
for such solutions affords a meaningful number of cases. 



5.2 Phase Space Boundaries 

Every point in the 3-dimensional P-V-9 space that results 
in a physical solution can be labeled by w^o of that solu- 
tion. The equation Wao{P,7 ,9) = C9, where C ^ 3 is some 
prescribed number, defines a 2-dimensional surface in the 
solution space. This surface is the locus of solutions whose 
velocity increases by factor C, i.e., every wind on this sur- 
face has Voo/vt ~ C. The volume enclosed by this surface 
corresponds to solutions with Voo/vt C, solutions outside 
this volume have Voc/vt < C. 

Figure ^ shows the projections onto the P-F plane of 
two such surfaces with representative values of Voo/vt- Each 
curve encloses the allowed region for winds whose Voo/vt 
is at least as large as the boundary mark. These bound- 
aries are best understood by considering the wind problem 
at a fixed F. Moving from the bottom of the figure along a 
line of constant F and increasing P, the curve for a given 
Voo/vt is first intersected at its lower branch. This intersec- 
tion defines a minimal P for physical solutions, which can be 
estimated by inserting 9 from the liftoff condition into equa- 
tion |62| a more accurate expression for this Pmin follows 
from the analytic solution for optically thin winds. The re- 
gion below the lower branch is forbidden for that particular 
Voa/vT- Increasing P further crosses into the allowed region, 
leading to winds whose Voo/vt exceeds the boundary mark. 
Eventually, quenching by reddening sets a maximum to the 
value of P that still allows Voo/vt as large as the boundary 
mark. This Pmax corresponds to the second intersection with 
the curve, and the region above the upper branch again is 
forbidden. The two branches of each Voo/vt boundary are 
defined by the two gravitational quenching modes and re- 
flect the central role of P in controlling both the dynam- 
ics and radiative transfer; this parameter sets the scales for 
both the acceleration (eq. and optical depth (eq. ps] ). 
The physical domain for P is between the two branches, 
Pmin < P < Pmax- As F decreases, the gravitational force 
increases and the two branches approach each other: Pmin 
becomes larger (it takes stronger acceleration to overcome 
gravity) and Pmax becomes smaller (it takes less reddening 
to quench the wind). The meeting point of the two branches 
defines an absolute minimum for F, determined exclusively 
by the grain properties. Below this minimum, liftoff requires 
such a high P that reddening quenches the wind immedi- 
ately. Note that this minimum is about 3-5 times larger 
than the Eddington limit F = 1, showing that the latter is 
not a sufficient condition for radiatively driven winds when 
the dust drift is taken into account (see below). 

The plots in figure |^ outline the proper phase space 
boundaries for radiatively driven winds, replacing the erro- 
neous bound Mvoa ^ L/c. If the velocity of a wind increases 
by a certain ratio, the parameters P and F of that wind nec- 
essarily lie inside the region bounded by the curve plotted 
for that ratio. The parameter 9 is not shown in this projec- 
tion but must conform to the bounds expressed by the liftoff 
condition and equation 
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Figure 7. Phase space boundary: The parameters P and F of 
winds whose velocities increase by the indicated ratios Voo/v^ 
must lie to the right of the correspondingly marked boundaries. 
The boundaries for amorphous carbon are marked by amC, for 
silicate grains by sil. 



5.2.1 The Lower Branch; Minimal M 

This branch reflects the constraints imposed by equation ^ 
and the liftoff condition. Equation ^ implies that Waa < P 
in all winds, namely t) < Up. As we discussed in ^ obser- 
vations usually give v <^ Vp. While this inequality poses a 
problem for the simple GS solution, the proper inclusion of 
drift and reddening shows that it actually provides impor- 
tant support for radiation pressure as the driving mechanism 
in dusty winds. 

From the liftoff condition (|2^), all winds must obey 
F > 1, the standard Eddington limit. Figure ^ shows that 
in dusty winds this bound is superseded by the combined 
effects of drift and reddening, which together impose a more 
stringent lower limit on F. The Eddington limit would be 
meaningful if the dust and gas were position coupled, in- 
stead the relevant regime in AGB winds is momentum cou- 
pling (see Oilman 1972). Indeed, F > 1 implies 



> 2.18x10" 

Mo 



(63) 



and this condition is obeyed by a large margin with typical 
parameters of observed winds. Although the Eddington limit 
does not impose a meaningful bound it provides another 
consistency check on the basic premises of the model. 
The liftoff condition also implies that 



> 



F + 3 
F - 1 



1/2 



(F-1) 



when F ^ 1 



(64) 

when F > 1 

Since F cannot approach unity (figure the meaningful 
constraint is 9 > F"'^. From its definition in equation GSl 
9 = M/Mc where 
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Mc = = 7.06x10' 



Q*L4 



1/2 



M 



yr 



(65) 



is a characteristic mass loss rate defined by the free param- 
eters. Therefore, the liftofi' condition implies a lower bound 
on the mass loss rate 



(66) 



When this constraint is violated, the grains are ejected with- 
out dragging the gas because the density is too low for effi- 
cient momentum transfer from the dust to the gas. 

The practical application of this result requires care. 
This constraint directly refiects our handling of the wind 
origin, the least understood part of the problem. Our liftoff 
condition is predicated on the assumption of positive ac- 
celeration right from the start. However, the initial accel- 
eration in fact is negative in the two-fluid formulation that 
neglects gas pressure; starting such a calculation with v = Vd 
produces a negative dv/dr, as is evident from both the NE 
(their equation 36) and HTT (eq. 7) studies. The wind is 
still lifted in these calculations because the dust accelerates, 
reversing the initial gas infall. But the negative derivative in 
turn reflects the shortcomings of these studies since the gas 
is in hydrostatic equilibrium in the absence of dust. This 
deficiency is removed in the study by Kwok (1975) which 
includes also gas pressure. At the sonic point Kwok finds 
TybaI T^-czM < 1 (his equation 15), in direct confiict with our 
liftoff assumption; the relation between the two forces is re- 
versed only after further grain growth. But Kwok's result, 
too, is unrealistic because it neglects the effect of radiation 
pressure on the gas molecules, likely to be instrumental in 
initiating the outfiow (e.g. Elitzur, Brown & Johnson 1989). 

An exact calculation of the minimal M is beyond the 
scope of our study. It would require proper inclusion of phys- 
ical processes, such as grain growth, whose understanding 
is still rudimentary at best. It is also instructive to recall 
that even the formulation of drag is only accurate to within 
". . . factors of order 2 or 3, which depend on accommodation 
coefficients and geometric factors" (Salpeter 1974). Fortu- 
nately, these uncertainties involve only the very initial stage 
of the outfiow and the severity of their impact is greatly re- 
duced once the fiow settles into the highly supersonic stage 
that is the main thrust of our study. Indeed, the minimal de- 
pendence on 9 of our final results demonstrates that for the 
most part, the initial conditions are quite irrelevant in the 
supersonic regime. While V^O > 1 remains a strict require- 
ment of our model, in reality this liftoff constraint could be 
less restrictive without an appreciable effect on most of our 
conclusions. 

We propose a phenomenological substitute for the 
mathematical liftoff constraint 



(67) 



where / (< 1) is an unknown factor. This only affects results 
that directly involve the liftoff relation — the location of the 
lower branch of the phase-space boundary (figure f?) and 
the estimate of the lower bound on M (equation p6). We 
estimate the phenomenological factor / by noting that it is 
mostly relevant at small P. In that region our solution gives 
(equation C4), which can be rewritten as 



r^e = f 



(68) 



With nominal values for vt and Vg, the requirement 
Voo/vt ^ 3 and consistency of the last two relations set 
/ ~ 0.1. As a result. 
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(69) 



replaces equation ^ as our phenomenological estimate of the 
minimal M. It may be noted that HTT proposed a similar 
relation. 



5.2.2 The Upper Branch; Maximal M 

The upper branch of the phase space boundary is defined by 
the quenching effect of reddening. Since it involves the wind 
outer regions where the supersonic flow is well established, 
it does not suffer from the liftoff uncertainties that afflict 
the lower branch. The upper branch deflnes a maximum M 
when the other parameters are fixed, discussed in paper II. 



6 GRAIN MIXTURES 

We now discuss the general case of a mixture of grains that 
can have different sizes and chemical compositions. For the 
most part, the results of the single-type case hold. A new 
feature is the variation of the mixture composition because 
of different drift velocities. 

The i-component of the mixture is deflned by its grain 
size ai and efficiency coefficients Q^.a; the latter may reflect 
differences in both size and chemical composition. As before, 
the details of production mechanism are ignored. We deflne 
the wind origin y — 1 as the point where the last grain 
type added to the mix enables the conditions for momentum 
coupling, initiating the supersonic gas outflow. We assume 
no further change in grain properties and denote the density 
of species i at that point ni{l) and its fractional abundance 
xii = ni(l)/nd(l). 



6.1 Scaling Formalism 

In analogy with the single-type case, the reddening profile 
of the i-th species is 

where 



^ J Qpr.ixBx{T,)dX. 

Since grain-grain collisions are negligible, different types of 
grains drift through the gas in response to the radiation 
pressure independent of each other. The velocity of species i 
is Vi and the ratio ni/nn varies in proportion to (,i{y) — v/vi. 
Introduce 



_9_ p4 

16^ 



qi = -p^ , where Q, = xi 



(71) 
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With the aid of the mixture average Q* we define Hm as be- 
fore (equation [ic|), and the drift function of the i-th species 



6 + y/wqi4>i 



(72) 



- qi4>i + y/wqi(j)i 
where w and 9 are the same as in the single-type case (eq. 



fad 



The radiation pressure force per unit volume is 



ni(y)naiq,(l)i{y). 



(73) 



We generalize the definition of cross-section area per gas 
particle (equation ^) via 

2 xu nd(l) 



/io-g = nai 



and 



/. = 1, 



(74) 



0(1) '^h(I) 

so that fi is the fractional contribution of species i to o-g. At 
every point in the shell runai — nno-gfid^i and the equation 
of motion is 



dy 



fiqi<l>i(y)Ci(y) - f 



(75) 



where P and F are defined as before. The equation of mo- 
tion retains its single-type form, the only modification is the 
weighted sum in the radiation pressure term. The dust prop- 
erties require now a larger amount of input but the other 
parameters remain the same — P, 6 and F. 

The radiative transfer equation requires as input the 
overall optical depth ry and dust density profile -q. Introduce 
the normalized density profile of the i-th species 





and the coefficients 



-^dy 

y^w 



Zt = 



^ Vi.fi I -J- dy 
/i y w 



(76) 



(77) 



where Qvi is the absorption efficiency at visual of the i-th 
species. Then 



(78) 



and the radiative transfer problem again does not require 
any additional input parameters. The most general dusty 
wind problem is fully prescribed by the grain properties and 
the three free parameters P, 9 and F. 

6.2 The Wind Velocity 

The wind velocity profile remains similar to the single-type 
case. Consider first optically thin winds, where — 1. An 
analytic solution does not exist now but the problem reverts 
to the single-grain case in two limits. When w is sufficiently 
small that w < qi for every i, drift dominates all grain types 
and the ^ = | profile of eq. BSl is recovered. In the opposite 
limit of w > qi for every i, drift is negligible for all grains 
and eq. ^ again applies with k — ^. In reality, the grain 
types are distributed between these two extremes. Different 
types move from the drift-dominated to the negligible-drift 
regimes as w increases, and the gas velocity profile evolves 
slowly from one shape to the other. Since the two profiles 



are quite similar, the overall behavior does not differ signif- 
icantly from the single-grain case. And since the only effect 
of reddening is to decrease k slightly further, the qualitative 
similarity remains for optically thick winds. 

6.3 Radial Variation of Grain Abundances 

Grains of different types have different velocities, creating 
a radial variation of the fractional abundances. The final 
abundance of species i is xnEi/E where 



Ci(oo) 
0(1) 



(79) 



and E = '^xnEi. Since ^ is monotonically increasing, 
Ei > 1 for every grain type and E > 1. The outflow en- 
hances the overall dust abundance, the final value of Ud/nu 
always exceeds its initial value. The fractional abundance is 
enhanced for grains with Ei > E and suppressed for those 
with Et < E. 

For simplicity, we discuss the behavior of this differen- 
tial enhancement only in the optically thin limit (0^ = 1), 
where Ei has the approximate behavior 



E, 



(I) 

\ 9 



1/2 



qi <9 

< qi < Wo 

< qi 



(80) 



This result is easy to understand. At a position where the gas 
velocity is w, grains with qt < w maintain a constant ratio 
Vi/v and have i^; = 1. Those with qi > w are falling behind 
and have = (w/qi)^^^. The expression for Ei reflects the 
two limiting cases of grains that are either moving with the 
gas {qi < 9) or lagging behind {qi > Woo) throughout the 
entire outflow, and the intermediate case of grains that are 
falling behind until the wind has accelerated to the point 
that w — qi. 

The differential enhancement is controlled by qi, the 
ratio of Q* for the particular grain type and for the whole 
mixture. The dependence of this quantity on grain size can 
be derived in a simple model that ignores spectral features, 
approximating the efficiency coefficient with 



Q{X,a) 



/'2naY 



A < 2na 



X > 2na 



(81) 



V A / 

This crude approximation, with ~ 1-2, is adequate for 
most purposes. Approximating the spectral averaging in the 
definition of Q« with its value at the peak wavelength of the 
Planck distribution produces 

{1 a > a, 

/„N/3 (82) 

where a, — 0.64 (i.m (1000 K/T*). It is interesting to note 
that r» = 2500 K gives a* = 0.25 [im, same as the upper 
bound of the MRN grain size distribution (Mathis, Rumpl 
& Nordsieck 1977). 

Consider an initial power-low grain size distribution of 
the form 



n^{a,y = 1) (X a ^ 



for 



^5 a < a* 



(83) 
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defined by the free parameters p and amin (the MRN dis- 
tribution corresponds to the special case with p — 3.5 and 
fflmin ~ 0.005 vim). With this general form, 



where ai — amin 



1//3 



(84) 



Only the lower bound amin enters here, q{a) is independent 
of the upper bound a*. The result for the abundance en- 
hancement factor is 



E{a) 



Vai6li/ 

{ e J 



0/2 
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a\Woo^^^ < a 



(85) 



The drift enhances preferentially the upper end of the size 
distribution, grains smaller than a certain size do not drift 
at all and their abundance is never enhanced. A power-law 
distribution that starts with index p emerges as a power-law 
distribution with the reduced index p — (5/2. To produce the 
MRN index of 3.5, the size distribution must start with a 
higher index of ~ 4-4.5. 

A complete theory for the grain distribution cannot be 
constructed without including the details of the production 
process. Our results show that some properties of the MRN 
distribution are natural features in the subsequent effect 
the outflow has on the fractional abundances. It is possi- 
ble that the MRN distribution could emerge as the only 
self-consistent solution of the complete theory. 



7 SUMMARY AND DISCUSSION 

Given grain properties, we have shown here that the su- 
personic phase of dusty winds is specified in terms of three 
independent dimensionless variables. A bank of solutions de- 
rived in the 3-dimensional space of these variables will con- 
tain the dimensionless density and velocity profiles of any 
possible wind around an evolved star. Matching the com- 
plete observations of any particular system requires at most 
scaling by an overall velocity scale. In addition, the solution 
of the radiative transfer problem is also contained in the 
outcome. We find that the solutions are generally controlled 
by a single parameter and can be described in terms of self- 
similarity functions that are independent of grain properties. 
Therefore, the results presented here contain the solution 
for the dynamics problem for all grains. The corresponding 
spectral energy distributions were presented elsewhere for 
amorphous carbon and silicate dust (IE95, IE97). Thanks 
to scaling, extending that database to any other grain com- 
position is as simple as a single run of DUSTY in which rv 
is varied over its entire feasible range. 

Our formulation starts once the dust properties are es- 
tablished. These properties, as well as the value of M, are 
set during an earlier phase of the outfiow. This initial phase 
is still poorly understood and our analysis has nothing to 
say about it. However, most observations involve the sub- 
sequent supersonic phase, which is the one addressed here. 
We associate the starting radius ri with prompt dust for- 
mation, but this does not enter directly into the solution 
because r\ drops out of our formulation. Extended, rather 



than prompt, dust formation can be accommodated with- 
out any change if ri is identified with the endpoint of grain 
growth. The association of r\ with dust condensation enters 
only indirectly — since the region r < ri is devoid of dust, 
the diffuse radiation vanishes at ri (IE97). Even this indirect 
effect is inconsequential in optically thin systems, where the 
dynamics impact of the diffuse radiation is negligible alto- 
gether. Therefore, our analytic solution eg. is valid in all 
optically thin winds irrespective of the nature of dust forma- 
tion; if grain growth is extended over a distance comparable 
with the dynamics length scale, the solution simply becomes 
applicable from the point where that growth stops. The only 
direct reference to the association of dust condensation with 
ri occurs when this radius is expressed in terms of other 
physical quantities (eq. |40| ). If the explicit dependence on 
ri is instead left intact, any arbitrary model of dust con- 
densation can be incorporated through the behavior it gives 
for r\. Young's correlation (sec. 4.2) implies that r\ oc L^^'^ , 



consistent with prompt dust formation. 

One of the grain parameters employed here is the cross- 
section area per hydrogen nucleus o"g. This quantity is di- 
rectly related to the ratio of dust-to-gas mass loss rates 



Md 4ps 



M 



3mn 



(86) 



where ps is the density of the dust solid material. This ratio 
(which requires additionally ps) was used as an input in the 
calculations of HTT, who denoted it 5 and pointed out its 
significance to the wind solution. Our analysis shows that (Tg 
actually does not enter into the scaling formulation of the 
problem; indeed, we have solved here the dusty wind prob- 
lem for two different types of grains without ever specifying 
the value of o-g for either of them. The only dust properties 
required are the efficiency coefficients Q\ and sublimation 
temperature Tc. Even the grain size is not strictly needed 
since it only enters indirectly as one factor, albeit a major 
one, in determining Qx- The dust abundance is an extrane- 
ous parameter controlling only the correspondence between 
the model results and physical quantities, just like the ve- 
locity scale which is irrelevant for the solution itself. 

Our analysis highlights the central roles that drift and 
reddening play in the velocity structure of dusty winds. In- 
dependent of grain properties, the significance of drift is 
diminished when P > a universal result refiecting the 
mathematical structure of the equation of motion. Redden- 
ing, on the other hand, becomes important at ry > 1 and 
the corresponding P is determined by Qy (equation [st]). If 
Qv were, say, 10~^ reddening would start playing a role only 
at P > 100, long after the drift ceased to be a factor. Since 
astronomical dust has Q\ of order unity, reddening starts 
controlling the force just as the drift dominance ends. 

Finally, why are typical velocities of AGB winds ~ 10- 
20 km s~^? Our results provide explanations for both the 
magnitude and the narrow range of v^o- While the driving 
force is radiation pressure, drift and reddening play a funda- 
mental role. Rapid drift at small ry decreases the dust abun- 
dance from its initial value, reducing the captured fraction 
of radiative energy. And at large optical depths this fraction 
is similarly reduced because the radiation spectral shape is 
shifted toward longer wavelengths where the absorption ef- 
ficiency is small. The two effects combine to produce the 
very narrow range of velocities displayed by equation ^ 
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— the range of variation with either L or rv is less than 
factor of 3, the dependence on dust properties is similarly 
weak. Given that the luminosity scale is ~ 10** Lq , the veloc- 
ity magnitude is determined by the dust parameters. Single 
grain properties cannot deviate much from the values listed 
in Table 1, the only parameter that is completely free is 
the dust abundance. Since the velocities are typically ~ 10 
km s~^, the dust abundance must be such that (T22 is of 
order unity; were that abundance 100 times higher, typical 
velocities would be ~ 100 km s~^ instead. This implies that 
the dust is produced with an abundance comparable to that 
found in the general interstellar medium and resolves the 
conflicts pointed out by Castor (1981). 
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APPENDIX A: GLOSSARY 



a 

Ei 



dust grain radius 

abundance enhancement of the i-th com- 
ponent of a mixture owing to differential 
drift velocities (eq. [zil) 



Ki...Ki 



L, Li 

M,Mo 

M,M-6 

rid, riH 

P = Vp/Vn 



ri,ri,i4 
T. 

Tc , Tc3 

7k, Tka 

V, Vd, Vrcl 
Vt 

Vac, VI 

W = v/Vm 
Woe = Vac /Vm 

y = r/n 
r 

5 = i/(r-i) 

C = v/Vd 





= Vt/Vu 



e 



(jg, 022 



TV 



similarity functions summarizing the so- 
lution observable implications (equations 
H, ^ and 1^) 

luminosity and its magnitude in 10^ L0 
stellar mass and its magnitude in M© 
mass loss rate and its magnitude in 10~® 
Mo yr-i 

number densities of dust grains and hy- 
drogen nuclei, respectively 
main scaling parameter of the dimension- 
less equation of motion; sets the ratio of 
radiation pressure to the drift effect 
the dust absorption efficiency at visual 
Planck average at the stellar temperature 
of the efficiency coefficient for radiation 
pressure (eq. ^ 

fractional contribution of the i-th compo- 
nent of a mixture to Q, (eq. ^l|) 
the wind inner radius and its magnitude 
in lO^"* cm 
stellar temperature 

the dust temperature at t\ and its magni- 
tude in 1,000 K 

the kinetic temperature at ri and its mag- 
nitude in 1,000 K 

velocities of the gas, dust and their differ- 
ence, respectively 

velocity scales characteri zing the radiation 
pressure (eq. |^, drift ( |lo| ) and gravity 
(^^, respectively 

the wind initial velocity (identified with 
the isothermal sound speed at ri; eq. ^) 
the wind final velocity and its magnitude 
in km s~^ 

dimensionless velocity 

the dimensionless final velocity 

dimensionless radial distance 

ratio of radiation pressure to gravity (eq. 

§) 

auxiliary quantity for the analytic solution 
of optically thin winds 
dimensionless drift profile (eq. ^|) 
dimensionless drift profile of the i-th com- 
ponent of a gas mixture (eq. ^2|) 
dimensionless dust density profile (eq. ^|) 
the dimensionless initial velocity 
the transformation function between P 
and TV (eq. ^) 

dust cross-section area per particle at con- 
densation and its magnitude in 10"'^^ cm^ 
(eq.|) 

overall optical depth at visual 
reddening profile (eq. ^) 

velocity- weighted harmonic average of 4> 
(eq. g 

density-weighted average of (eq. |42| ) 
a dimensionless function determined by 
radiative transfer and its value in optically 
thin winds (eqs. ^ and ^l|) 
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APPENDIX B: DRIFT VELOCITY 

The drag force on a grain moving at velocity Wrei through 
gas with density p is 



-Fdrag = vra p X 



L V: 



1 when v^ci > 
ivt when Vici < 



(Bl) 



for the subsonic {vici < vt) and supersonic (wrci > vt) 
regimes (e.g., Kwok 1975). The drift reaches steady state 
when -Fdrag ~ fpr, the radiation pressure force on the grain. 
Introduce 



,,2 _ -Ppr 



TiaP-p 



cp 



(B2) 



then the steady-state drift velocity is 
vu when vu > vt 

^ when vo < vt 



Vrcl 



(B3) 



Both limits can be combined in the simple expression 



Wrcl 



Vd + VT 

Kwok proposed instead the more complex expression 



{l [(44 



2 

Vt 



1/2 



(B4) 



(B5) 



He derived this result by combining first into a single ex- 
pression the two limiting forms of the drag force, instead of 
Vj-ci itself, and then solved the resulting quadratic equation 
for VtbI- Kwok's expression and ours are based on the same 
ingredients and give identical results in the limits of both 
Vd > Vt and vu < vt- 



APPENDIX C: ANALYTIC APPROXIMATIONS 



This expression is an expansion in the problem's small pa- 
rameters. Consider first the limit (3 = 0, which gives 



4 3/2 
— ?/? 

3 "^oo 



(C3) 



At small values of Woc the first term on the right dominates, 
at large values the second. Keeping only the dominant term 
produces = {j^P* f^^ (1 - e3/2-)-2/3 -^^ ^^le first regime, 
Woo = ^'(1 — e^)~^^^ in the second. Neglecting the e cor- 
rections, the two terms and the approximate solutions they 
yield are equal to each other at ~ P ~ ^ ■ Since small 
uioo ensures ui < 1 for the entire wind, C, = w^/^ in that 
case whileC = 1 for most of the wind when w^o 2> 1 (see 
equation |24[ ). Therefore, P = ^ is the transition between 
drift dominance and negligible drift, and 



(^P*) ' for P < 16/9 (drift dominated) 
P for P > 16/9 (negligible drift) 



(C4) 



These two limit expressions for Wod can be combined in the 
single form listed in equation ^ 

Finite (3 corrections are generally small when e ^ i 
and P ^ ^. The deviations from unity of the terms inside 
the large square brackets are less then 30% in that domain. 
The only significant corrections can come from the terms 
in /Je^/^wi^^ Since < P, Pe'-^'^wU,^ < 0.3P^^^ for aU 
e ^ I and /3 ^ i. Therefore, this term is always negligible 
in the drift dominated regime. But at large P it can become 
significant, at P = 10 it already amounts to a 90% cor- 
rection. When this term dominates, the solution becomes 
u)^ = P'^/0^e. However, reddening effects become signifi- 
cant before this limit is reached, therefore we are justified in 
maintaining the /? = approximation; figure |l] shows that 
this faithfully presents the entire e =^ | and (3 ^ \ domain. 

We proceed now to the velocity profile. With /3 = 0, 
equation M gives 



The quantities of interest obtained from the solution of the 
equation of motion are the profile shapes of the velocity and 
density variation, the final velocity Woo and the optical depth 
ry. Introduce 

5 



f3- 



(Cl) 



The profile shape is conveniently expressed in terms of u, 
which varies from e to 1, with e (< i) and [3 (< 1) as small 
parameters. Furthermore, we note that 9, too, is a small 
parameter in all practical cases and use 5max = 9 + \fQ — \fQ 
as well as (5 — 6 ~ /3v^. 



CI Negligible Reddening 

When 0=1, the complete solution (eq. ^3[) gives the fol- 
lowing equation for Wao'- 



-l/2,„l/2\ 



3fl 1/2 



l + f/3. 



3/3^6 -£^''^(1 + 1/? + 3/3^) 



,3 3/2, l-Pe^'^ 



+ 3r.-ln^_^^ 
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(C2) 



— + 



(^ 



.3/2^ 



1 - e2 + %w[ 



-1/2 



(1 



.3/2-) 



(C5) 



Similar to equation OS, this equation changes its behavior at 

-1/2 



; at smaller Woo the terms proportional to 



dominate, at larger Waa the other terms. The two limits give 

r 2/3 for P < 16/9 
k = \ (C6) 
il/2 for P> 16/9 



u=\l- 



i/fc 



y 



From this we can immediately determine the dust density 
normalized profile rj (oc (/y'^u), required for the radiative 
transfer equation. The ratio (^/u is proportional to 
in the drift-dominated regime and to 1/m when the drift is 
negligible. With the particular value for k in each case, (/u 
is the same as 1/u with the power index —k replaced by 
1 — k, leading to 



A 

y 



y 



y 



A 



i/k 



1 



(C7) 



The velocity and density profiles are independent of P, 
which is relevant only for the choice of the solution regime. 
The only parameter that enters explicitly is e, and even 
this dependence is confined mostly to the wind origin where 
u(l) = e and 77(1) = A/e'^~'°'/*. However, the e-dependence. 
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which must be maintained to avoid a singularity for rj at y 
= 1, rapidly disappears once y > 1 + e. 



C2 Reddening Corrections 

Reddening corrections are conveniently expressed in terms 
of the quantity <1>. Since <j!> decreases away from the wind 
origin, the averaging is dominated by its upper end and 



$ ! 



<t>{y 



1 



(C8) 



For simplicity, only 5 = is considered. The equivalent of 
equation C3 is then 



4 



3/2 



du 

7f 



^ du 



(C9) 



Approximating the integrand in the first term by its largest 
value, our approximate solution for Woo becomes 



A 3/2 2 
3%/¥ * 



In analogy with equation Cc this yields 



for P < f$ 



16^1/2 



for P > 



(CIO) 



(CU) 



C3 Optical Depth 

The relation between P and rv is obtained by combining 
equations BS^ and IClCl. When (5 = 0, 



Qv 



3/2 



+ 



TV 

Qv 



(C12) 



From this we can find the value of P that yields ry 
With the approximation $ ~ 1 at this point we get 



P = 



Qv V 



1 + 



.1/2 



1/2 



(C13) 



Finally, the (5-dependent correction term in equation ^ can 
be estimated in the optically thin case. With the leading- 
order velocity profile from equation |C6| , the integration re- 
sult is 



dy 



1 1 



^(i-k)/k 
- ei/fc • 



(C14) 



<p(— 
y 



(D2) 



Given (j), this equation is solved by iterations that start with 
the analytic solution for w derived in the previous section. 

The code DUSTY (Ivezic et al, 1999) starts by solving 
the radiative transfer equation with an initial density profile 
»7 taken from equat ion bTl . The resulting reddening profile (f) 
is used in equation [D1| to find w. A new reddening profile is 
calculated with the corresponding r) (see equation E5|), and 
the process is repeated until (j>, V and w are self-consistent. 
Convergence is rapid and the number of radial grid points 
is modest — less than 30 points are adequate for optical 
depths up to 100. The power of the method stems from the 
fact that the steepness of the velocity and density profiles is 
built in right from the start. 

We find it advantageous to replace P, F and 6 with 
three other independent input parameters. First, because of 
the central role of radiative transfer we specify the optical 
depth TV and determine P from 



TV 

NQv 



where A'^ — 



c 



dy 



(D3) 



(D4) 



(see equation |25|). Next, consider the quantity 
1-i 

_ y 
z{y) 

the maximum of the prescribed profile anywhere in the shell. 
The wind stalls when (/max ~ F, therefore physical solutions 
have (7max ~ /F with / < 1. We choose / as an input pa- 
rameter instead of F with the replacement 

1 _ / 

r C'max 



(D5) 



The parameter / varies from for negligible gravity to 1 
for gravitationally quenched winds. Finally, instead of 9 we 
specify as an input parameter e = 0/wca- If Y denotes the 
shell outer radius, 9 is derived from the formal solution at 
that point, which gives 



z{Y) 



(D6) 



The advantage of tv, e and / over P, 9 and F as input 
parameters is the more direct contact they have with phys- 
ical quantities. In particular, the choice e ^ | and / < 1 
guarantees a priori a physically meaningful solution. 



APPENDIX D: NUMERICAL PROCEDURES 

The velocity and density profiles are very steep near the 
wind origin (cf eq. |C7| ). As a result, the solution of the equa- 
tion of motion e4 can become a difficult numerical problem 
that requires a prohibitive number of radial grid points, es- 
pecially when optical depths are large. To bypass these prob- 
lems, we find the velocity profile instead from the integral 
equation expressed by the formal solution 



^9^ + 



(Dl) 



where 
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